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£h Abstract 

We study fractional exclusion statistics for quantum systems with SU(2) Byrn- 
es ' 

metry (arbitrary spin S), by generalizing the thermodynamic equations with 

•rH . 

squeezed strings proposed by Ha and Haldane. The bare hole distributions 
as well as the statistical interaction defined by an infinite-dimensional matrix 
specify the universality class. It is shown that the system is described by the 
\evel-2S WZW model and has a close relationship to non-abelian fractional 
quantum Hall states. As a low-energy effective theory, the sector of mass- 
less Tiis parafermions is extracted, whose statistical interaction is given by a 
finite-dimensional matrix. 
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Fractional exclusion statistics proposed by Haldane ]L) has attracted much current inter- 
est. In particular, the notion of ideal exclusion statistics is useful to characterize elementary 
excitations, which naturally interpolates free fermions and free bosons Jl]-§J. The ideal 
exclusion statistics is featured by statistical interaction g^ v defined as, 

dD„ 



ON,, 



-g»v, (i) 



where g^ u depends on the species /i and z/, but not on the momentum kj. Here and iV M 
denote, respectively, the number of holes and particles with species /i. A typical system in 
which ideal exclusion statistics is realized is quantum 1/r 2 models in one dimension |)]||. 

The ideal exclusion statistics has a close relationship |7|,[8| to abelian Chern-Simons the- 
ory for the fractional quantum Hall effect (FQHE) ||. It turns out that the statistical 
interaction matrix plays the same role as the topological-order matrix in FQHE, and hence 
the multicomponent model in ref. may serve as a model for the edge states of the FQHE, 
namely, its low-energy behavior is described by multiple c = 1 conformal field theories 
(CFTs). In this sense, we may call such statistics abelian fractional exclusion statistics. In 
theories of the FQHE, another interesting class of FQH states has been proposed, in which 
quasiparticles obey non-abelian statistics |TI]-[rjJ. The corresponding theory is given by, 
for example, non-abelian Chern-Simons theory, and is closely related to c > 1 CFT. From 
these observations, a natural question arises, i.e. how can one generalize fractional exclusion 
statistics to non-abelian cases? 

We address this question in this paper, by applying the idea of exclusion statistics to 
SU(2) higher-spin models in one dimension. We propose basic equations for ideal statistics 
with arbitrary spin based on the idea of squeezed strings introduced in the S = 1/2 model 
T5| . We then study the thermodynamics and show that the critical behavior is described by 



the level-25 1 SU(2) WZW model. It is found that in order to formulate non-abelian statistics, 
the bare hole distribution also plays a vital role in addition to the statistical interaction. 
We next deduce explicitly the statistical interaction of massless Z2S parafermion which 
characterizes the nontrivial properties of non-abelian statistics. The relation to non-abelian 



FQHE is mentioned briefly. 

To begin with, let us briefly review the ideal spinon description of the quantum system 
with S = 1/2 [|TJ. In this case, ideal statistics is modeled by lattice bosons with statistical 
interaction g = 2, reflecting level-1 SU(2) Kac-Moody symmetry. Recall that the ideal 
statistical interaction g^ u is simply related to the two-body phase shift of the scattering 
matrix in one dimension through — kj) = n (g^ — 5^ u )sgn(ki — kj) fjjfl,!]]- Note that the 
phase shift with step function is inherent in ideal statistics. Thus the two-body phase shift 
in the S = 1/2 model reads <fi = 7r sgn(A;j — kj). Using this phase shift one could naively write 
down the Bethe equation to determine the spectrum for the periodic chain. The resulting 
equation, however, does not yield the right number of eigenstates. To overcome this difficulty 
Ha and Haldane introduced the fictitious string solutions called "squeezed strings" fc^eR 
(/i = 1,2, ••■) [p~5| . The Bethe equation now becomes the coupled equations for strings. 
What is remarkable is that these string solutions give not only complete degeneracy of 
states but also the correct eigenvalues even for finite systems. This method thus provides 
precise thermodynamics of the free spinon gas obeying ideal exclusion statistics. 

We now wish to generalize this description to the higher spin model. We expect that 
the basic equation for ideal statistics may have a similar structure to the Bethe equation for 
integrable models. Examining carefully the Bethe equation for solvable Heisenberg chain 
with higher spin [|]1||17]] , we are led to propose the following basic equation, which generalizes 
the Ha-Haldane equation to arbitrary spin-,!?, 

d^Nk? = 2-kI? + vr ]T <>sgn(tf - fcj) (2) 

for /i = 1, 2, • • •, where d>~> = min(/i, 2S), c^ u = G^ u — 5^ u with 

G^v = 2min(/i, u), (3) 

and N is the number of sites. Each k is defined in the region k e [— tt, 7r]. In the above 
equation, we have introduced the squeezed string solutions kj, which are coupled with each 
other by the matrix c MI/ . The statistical interaction of this model is given by the infinite 
dimensional matrix G. The total energy is given by E = Yl/jt,i e ° (^f)j where 
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el(k) = d^e°(k), m = \{e-^). (4) 

Here we have assumed the quadratic dispersion of the energy, which is indeed the case for 
integrable 1/r 2 models ||, though the following analysis is not sensitive to the form of e°. 
Note that the r.h.s. of eq.(^) takes the same form as the Ha-Haldane equation. Thus the 
effect of higher spin is incorporated through the factor S s > in eqs.(|2]) and (|J), whose role is 
crucial to identify the universality class of the system. We will show momentarily that 
is related to the bare hole distribution. 

In the thermodynamic limit, eq.(H) leads to the coupled equations for the distribution 
functions of particles (p) and holes (p^), 

PM + pf{k) = J-dj?> -£c^(fc). (5) 

The equilibrium distributions are then determined by minimizing the free energy SF = 0, 
which gives 

e^k) = eJ(A;) + pH + T £ c, u ln(l + e"-^), (6) 

V 

where H is a magnetic field, and the dressed energy e is introduced by p^ h '/p = e £//T . The 
free energy per site / = F/N at equilibrium is then given by 

f = -HS-^-Y. d i S) r dkln(l + e-^ T ), (7) 

ZTT " J-n 

A 1 

which can be rewritten by using eq.(^), 

Eqs.fl^), (|6|) and (||) are the basic equations for describing thermodynamics of the systems 
obeying non-abelian exclusion statistics with an arbitrary spin S. 

Let us now examine low-temperature properties. First, by using G~l = (25^ — 6^ u+ i — 
<Vfi,i/)/2, ec L-(H) can be converted into 

el e° 



where L M = ln(l + e €M//T ). We can see from this equation that > 0, while e 2 s can change 
its sign at some value k = ±Q (pseudo- Fermi points). At T = 0, eq.(|) reduces to 



e li (k) = ^(k)+ t jiH-G ia se^\k) for M ^ 2S, 
4s (*0 = e° 2S (k) + 2SH - G 2S2S e i 2S ) (k), 



(10) 
(11) 



where = 6(e 2S (k)) and e$ = -9(-e 2S (k)), with 9(x) = 1 (0) for x > (x < 0). 

We start with the case of H = to simplify the discussions. At T = 0, we immediately 
find 



S, 



(12) 



which implies that the ground state consists of the sea of 25'-strings with the pseudo- Fermi 
points ±Q with Q = n. At low temperatures, we can thus see from (§) that the most 
important contributions to the free energy come from the vicinity oik = ±Q. For example, in 
the region close to the right Fermi point, [it — e, tc], the driving term in eq.(0) is approximated 
by ir(k — n)/2T = —e~ d where 9 is defined in the region 9 G [— ln(7re/2T), oo). In the limit 
T — > 0, we can set 9 G (—oo, oo). We get from eq.(^) 
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(13) 



where (p^(9) = e^/T. The free energy is then given by 

£ / d9(^L, - ^l; 

1 J —oo 



4vr 2 



(14) 



where we dropped the contribution from the ground-state energy. In order to evaluate this 



integral, we must solve eq.(|T3"D in the asymptotic region 9 — > ±oo. We find |L7 

for 1 < /I < 2S - 1, 



{ ( . (/i + l)7T\ 

sin - 



2S + 2 



IT 



sin 



V 25 + 2 J 
(fi-2S + l) 2 
(/i+1) 2 , 



for 25* < /i, 



(15) 



(16) 



where = e^^ 00 ^ + 1. Consequently we have the low-energy expansion of the free energy, 

T 2 25 ttT 2 3k 

where the factor 2 in the first expression is due to the contributions from both the left and 
right Fermi points, C(x) is Rogers dilogarithm function defined by 



ln(l -*) Int 



t 1 - t 



(18) 



k = 2S and v = it is the Fermi velocity of the 2S'-string sector. Here we have used dilogarithm 
sum rules [fLS ]. According to the finite-size scaling law |T9|j20[| , we see that the critical 



behavior of the present system can be described by the SU(2) 2 s WZW model. Namely, 
statistical interaction G in eq.(|) with additional in eqs.(@) and (|J) results in the central 
charge, 

3k , ~ 

C -—2- (19) 

with k = 2S. This result implies that the universality class of the present system is charac- 
terized not only by the statistical interaction G but also by the factor S s \ Hence, also 
plays a role in formulating ideal statistics for general S. Let us check how enters in the 
definition of exclusion statistics ([I]). To see this, we integrate ([I]) to obtain [|| 

D^ = -Y,9^N U + Dl. (20) 

V 

One can see that is directly related to the integration constant _D° such that = 
Nd^/2n. This term can be interpreted as the bare distribution of holes when there are 
no particles in the system. Therefore, the bare hole distributions as well as the statistical 
interaction are the fundamental elements which characterize fractional exclusion statistics 
for non-abelian cases. This point has not been seriously taken into account so far, because 
the above problem does not happen for abelian cases. 

We now derive low-energy effective theories. In the low-energy region, we can explicitly 
separate the theory into the abelian c = 1 sector and the c = 2(k — l)/(k + 2) sector. 
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This analysis enables us to naturally formulate ideal statistics with Z25 symmetry. We will 



extract the Z 2 s parafermion sector following Tsvelik pT|. First, we switch on magnetic fields 



to deduce effective equations for modes 1 < /i < 2S — 1. At T = 0, we have from eq.flTO), 



and from eq.flll 



TTT^as for »<2S, 

e,(k) = {2S* S (21) 

(/1 - 2S)H + eY s > for 2S < //, 



, 2S(e°(k) +H) for Q <\k\< n, 

e^{k) = { M - (22) 

for < Q, 

■ for Q < \k\ < n, 

eis\k) = { 1 (23) 

-(e°(k)+H) for \k\<Q, 



where e 2 s(^Q) = with Q = \lix 2 — 2H. From these equations we see that the p, < 2S 
modes (parafermions) as well as fi = 2S mode (U(l)) are still massless, while those for 
2S < fi get massive. This is a novel feature which has not been observed for the integrable 
spin chain |TB|,|TT| in which all modes become massive under magnetic fields except for 



2S"-string sector |HJ]. The existence of massless parafermions is closely related to ideal 
statistics, i.e. it directly reflects the step-function of the phase shift. If the deviation from 
ideal statistics is introduced, the parafermion sector immediately becomes massive, as will 
be discussed later. In this sense, the present model is regarded as an ideal model for the 
massless spin chain with arbitrary S. 

Let us proceed to calculate the distribution functions. Subtracting eq.@ for fi < 2S 
from that for fi = 2S, we obtain 

25-1 



with ji, v = 1, 2S — 1, where 



$ + Pm = §P^5 j - E (GW - «fr ( 24 ) 



Gfj,v — — — . (25) 
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The modes \i < 2S still couple to the 25'-string sector. Since the U(l) and Z 2 5 sectors are 
both massless, it seems not easy to separate them. Nevertheless, if we replace p 2 s i n the 
above equation by the ground-state distribution, 

S/ir for Q < \k\ < tt, 
PS = II—' (26) 

I for \k\ < Q, 

we then get desirable equations describing Z 2 s parafermions from eq.(|24]) with the dis- 
persion fl2~T|). This can be easily confirmed by examining the thermodynamics again 
based on these equations and computing the central charge, which in fact results in 
c = 2{k — l)/(k + 2). Hence we conclude that fl2"l| ) with fl2"E| ) describes ideal exclusion 
statistics for Z 2 s parafermions. Their statistical interaction is given by (2S — 1) x (2S* — 1) 
matrix G in eq . (f25|) . Note that G is the inverse of the A 2 s Cartan matrix C: 



CG = 21. (27) 

Therefore it can be expressed by the incidence matrix I; G — I = 1(2 — Z) -1 . Having this 
expression here is natural since it has already appeared in the thermodynamic Bethe ansatz 
for the massive perturbation theory of CFT whose ultra-violet behavior is described by Z N 



parafermions [22 



It is instructive to consider a question why the present system has massless parafermions 
even in magnetic fields in contrast to the ordinary solvable spin chain models [|T^,|T^,^] . To 
clarify this point, let us take equations, 

dl S) Nk(\?) = 2nl? + ~ A;), (28) 

where k(X) = 2arctanA, = 0\^ u \ + 29\ IM _ u \ +2 + ... + 2# A1+J ,_ 2 + 6^ +v (for fi = v the first 
term is omitted) with ^(A) = 2 arctan(2A/r/i). This equation recovers (§) in the r — > 
limit. Also, one notices that r.h.s. is the same as the Bethe equation for the solvable XXX 
spin chain for T = 1. Using this equation, we can derive dressed energies for fi < 2S modes, 

e,(X)= f°° dX'K.iX-Xy^iX'), (29) 

J —00 



where 

K (x] _J smjnfx/N) 

M ' AT cosh(27rA/AT) + ca&^n/N) ' 1 ' 

apart from a term which disappears with r — ► 0. For T = 1, this expression coin- 
cides with that for the XXX spin chain derived by Tsvelik, giving a massive dispersion 
e At (A) oc sin(7r/i/iV) cosh(27rA/A^) for small A |[21|| . Contrary to this, in the r — > limit, 
K^{\) — > (fj,/N)S(X), and this gives massless dispersions fl2lD . One can thus see that mass- 
less parafermions are realized only for ideal statistics, and infinitesimal deviation from it 
may bring about massive parafermions in magnetic fields. 

Finally, let us briefly mention the relation to the edge states of the FQHE. As proposed 
by Moore and Read fL0| , and Block and Wen [fTT|-|l4}|, there are several models for non- 



abelian FQH states. An example is the FQH state with U(l)x SU(m) n symmetry fL2,14 
The present results can describe the behavior of the non-abelian SU(2) n sector of the above 
edge states. Another interesting example is the FQH state with XJ(l)xZ N model proposed 



by Moore and Read flit) JIT , D5|,H|1. The "Pfaflian" part of Zn is essential for this non-abelian 
FQH state. Our result for Z 2 s symmetry extracted from the SU(2) 2 ,s case may describe the 
above sector in the edge state. It may be interesting to generalize the present discussions to 
describe non-abelian statistics with more complicated symmetries. 

This work is supported by the Grant-in- Aid from the Ministry of Education, Science and 
Culture, Japan. 



9 



REFERENCES 

JSPS Research Fellow. 
[1] F. D. M. Haldane, Phys. Rev. Lett. 67, 937 (1991). 

[2] M.V.N. Murthy and R. Shankar, Phys. Rev. Lett. 72, 3629 (1994); preprint IMSc-94/24 
|cond-mat/94U4U96| . 

[3] Y.-S. Wu, Phys. Rev. Lett. 73, 922 (1994). 

[4] D. Bernard and Y.-S. Wu, preprint SPhT-94-043, UU-HEP/94-03, |cond-mat/ 9404025 . 



[5] F. Calogero, J. Math. Phys. 10, 2197 (1969): B. Sutherland, J. Math. Phys. 12, 246, 
256 (1971); Phys. Rev. A4, 2019 (1971); A5, 1372 (1971): F.D.M. Haldane, Phys. Rev. 
Lett. 60, 635 (1988): B.S. Shastry, Phys, Rev. Lett. 60, 639 (1988): Y. Kuramoto and 
H. Yokoyama, Phys. Rev. Lett. 67, 1338 (1991). 

[6] See, for reviews, N. Kawakami, Prog. Theor. Phys. 91, (1994) 189; 

F. D. M. Haldane, in Correlation Effects in Low Dimensional Electron Systems, eds. A. 
Okiji and N. Kawakami (Springer, Berlin, 1994). 

[7] T. Fukui and N. Kawakami, Phys. Rev. B51, 5239 (1995). 

[8] T. Fukui and N. Kawakami, preprint YITP/K-1107, |cond-mat/9504118 

[9] X.G. Wen and A. Zee, Phys.Rev. B46, 2299 (1992). 
[10] G. Moore and N. Read, Nucl. Phys. B360, 362 (1991). 
[11] B. Block and X.G. Wen, Nucl. Phys. B374, 615 (1992) 
[12] X.G. Wen, Int. J. Mod. Phys. B6, 1711 (1992). 
[13] X.G. Wen, Phys. Rev. Lett. 70, 355 (1993). 
[14] X.G. Wen, preprint MIT-CMT/June 95, |cond-mat/9506066 



10 



[15] Z.N.C. Ha and F.D.M. Haldane, Phys. Rev. B47, 12459 (1993). 
[16] L.A. Takhtajian, Phys. Lett. 87A, 479 (1982). 

[17] H.M. Babujian, Phys. Lett. 90A, 479 (1982); Nucl. Phys. B215, 317 (1983). 



[18] A.N. Kirillov, preprint , |hep-th / 9408 1 1 3| : L. Levin, Dilogarithm and associated functions, 



North Holland, Amsterdam, 1981. 
[19] H.W. Bldte, J.L. Cardy and M.P Nightingale, Phys. Rev. Lett. 56, 742 (1986). 
[20] I. Affleck, Phys. Rev. Lett. 56, 746 (1986). 

[21] A.M. Tsvelik, Nucl. Phys. B305, 675 (1988); Sov. Phys. JETP 68, 1436 (1988). 
[22] T.R. Klassen and E. Melzer, Nucl. Phys. B338, 485(1990). 



11 



